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Abstract 

We present results of extensive statistical and bit level tests on three imple- 
mentations of a pseudorandom number generator algorithm using the lagged 
Fibonacci method with an occasional addition of an extra bit. First imple- 
mentation is the RCARRY generator of James, which uses subtraction. The 
second is a modified version of it, where a suggested error present in the 
original implementation has been corrected. The third is our modification 
of RCARRY such that it utilizes addition of the carry bit. Our results show 
that there are no significant differences between the performance of these three 
generators. 
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1 Introduction 



Random numbers are needed in various applications, including cryptography M, 
stochastic optimization jjj], and Monte Carlo methods [|J. Because of practical 
reasons random numbers are usually produced by deterministic rules, implemented 
as pseudorandom number generators. In spite of their fully deterministic origin the 
quality of pseudorandom numbers may often be good enough for many applications. 

To confirm the suitability of a given pseudorandom number generator for practical 
use, it should be subjected to a rigourous test program which reveals the strengths 
and weaknesses of the algorithm and, in particular, its implementation. Recently, 
such an extensive test program has been carried out by the present authors jni . By 



performing a comparative evaluation using statistical, bit level and visual tests we 
were able to assess the quality of a group of random number generators which are 
commonly used in the applications of physics. 

One of the generators included in Ref. [P| was RCARRY, which uses the so called 
' 'subtract- and-borrow" algorithm which has been implemented by James ||. In the 
tests, RCARRY clearly displayed the poorest statistical properties of the genera- 
tors tested, suggesting possible problems in the implementation. Supporting this, 
James has recently reported the observation of M. Luscher that the original im- 
plementation of RCARRY may contain a small error, which may adversely affect 
the quality of the random number sequence. The purpose of the present work is to 
address this issue. To this end, we present results of extensive statistical and bit 
level tests on the corrected version of RCARRY, and compare the results to those of 
Ref. In addition, we test a slightly different version of the RCARRY generator, 
which uses an "add- and- carry" algorithm based on the addition of a carry bit. We 
call this generator ADCARRY. Our results reveal that there is very little difference 
between the statistical properties of the original RCARRY and its corrected version, 
as well as the ADCARRY generator. All these generators display a relatively poor 
performance in two of the gap tests presented here. 



2 Implementation of the Generators 

The three pseudorandom number generators tested in this work are based on a lagged 
Fibonacci algorithm, which is augmented by an occasional addition of a carry bit. 
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The basic formula is: 



X % = (X,_ 24 ± X t _ w ± c) mod b. (1) 

The carry bit c is zero if the sum is less than or equal to b, and otherwise "c = 1 
in the least significant bit position" ||. The choice for b is 2 24 . The period of the 
generator is about 2 1407 || and it produces random numbers distributed between 
[0,1). Only the 24 most significant bits are guaranteed to be good. 

The inclusion of the carry bit c in the lagged Fibonacci algorithm was done in order 
to improve its properties [16]. Recently, however, it has been shown [||, [L7j that this 



type of algorithms are in fact equivalent to linear congruential generators with very 
large prime moduli. Consequently, they inherit unfavourable lattice structures in 
higher dimensions. 

The original implementation of Eq. (1) was done by James ||, based on the ideas 
of Marsaglia et al. [0. It uses the subtraction contained in Eq. (1). In this work, 
we shall denote it by II. The second generator 12 includes the suggested correction 
of Luscher and James, who recommend replacing line 13 of the code of Ref. || 

uni = seeds(i24) - seeds(j24) - carry, 

by 

uni = seeds (j 24) - seeds (i24) - carry. 

The third generator ADCARRY (13) uses the operation known as "add-and-carry" , 
in which subtraction in Eq. (1) has been replaced by addition. In this version the 
lines 13 - 15 of are rewritten as: 

uni = seeds (j 24) + seeds (i24) + carry 
if (uni .ge . 1 . ) then 

uni = uni - 1 . 



Otherwise, the implementation is identical to that of RCARRY flT2", [13] 



3 Test methods 



Tests scrutinizing the quality of random numbers can be divided into three main 



categories: statistical tests |I| , bit level tests 
random numbers on binary level, and visual tests 



15| , |l"8H for testing the properties of 
| which may give some further 



qualitative information on the statistical properties of random numbers. A number 
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of these tests were implemented and employed extensively in Ref. ||18|| . In this work, 
we have repeated the same statistical tests for 12 and 13. They are listed in Table 
1, where the numbering refers to the parameters of Ref. |Tj|. From bit level tests, 
only the <i-tuple test [[|, [15] was done since it was shown to be sufficient. Finally, the 
random numbers were plotted in two dimensions for purposes of visual inspection. 

The test bench is described in detail in Ref. ]nj. Description of the statistical tests 
can also be found in Ref. [ 10| . In brief, the statistical accuracy of all the tests was 



improved by utilizing a one way Kolmogorov - Smirnov test |L0| to a large number 
(1000 or more) of test statistics. This approach has been realized earlier by L'Ecuyer 
[ |llTl . The final test variables are therefore the values K + and K~ of a Kolmogorov 
- Smirnov test statistic K |10[ . In each test the generator was considered to fail the 
test if the observed descriptive level 5 = P(K < {K + , K~}\H ) was less than 0.05 
or larger than 0.95. 



4 Results 

Results of the statistical tests for the descriptive levels 5 + and 5~ are summarized 
in Table 2, where the numbering refers to Table 1. In each test the chosen generator 
was initialized with the seed 667790. In case a failure occurred, the generator was 
subjected to another test starting from the final state of the first test. If another 
immediate failure occurred, the generator was tested for the third time starting from 
a new state with an initial seed 14159 (from the decimals of 7r). 

In Table 2, frames with thin lines indicate a single failure, frames with double single 
lines two consequtive failures, and frames with bold lines three consequtive failures 
in the corresponding tests. The results of the original RCARRY (II) by James || 
are shown on the left (from Ref. [|18|]), whereas the results of the corrected version 
(12) and ADCARRY (13) are at the center and on the right, respectively. 

Based on the results, it is clear that the corrected version of RCARRY (12) using 
arithmetic subtraction performs no significantly better than the original RCARRY 
(II). The main malady of RCARRY, namely the clear failing of the gap tests 6 and 
8 with parameters a — 0, j3 — 0.05 and a = 0.95, /3 = 1 HlOl , respectively, is still 
characteristic of 12. This signifies the existence of local correlations in the vicinity 
of zero and one. The same conclusion applies to ADCARRY as well, signaling basic 
problems with these algorithms. 
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In the (i-tuple test each implementation was tested two times and the bits considered 
failed had two consequtive failures. The results are shown in Table 3. In our 
notation, bit number one is the most significant bit (excluding the sign bit). For 
the original implementation of RCARRY (II) only the 24 most significant bits are 
guaranteed to be good, which the tests confirm |Tj|. The implementation 12 yields 
identical results, whereas ADCARRY (13) gives only 22 good bits (see Fig. 3). 

Finally, visual tests on bit level support the results above. In Figs. 1, 2 and 3 we 
show subsequent random numbers for II, 12 and 13 in binary form on a 120 x 120 
matrix, when only 24 most significant bits are included. No clear correlations are 
visible, except for the last two bits of ADCARRY where strong correlations are 
apparent. No visual indications of the suggested |5], [I7| lattice structure were found 
in these generators. 



5 Summary and conclusions 

In this work, we have compared the results of detailed statistical and bit level tests 
for three implementations of random number algorithms using a lagged Fibonacci 
sum with the addition of a carry bit. Results for RCARRY and its corrected version 
show very little difference. Also, a new generator ADCARRY using purely additive 
arithmetics fares no better statistically, and has two good bits less than RCARRY. 
Fortunately enough, these bits are not among the most significant ones. Overall, 
our results suggest that the basic algorithm of Eq. (1) on which these generators are 
based seems to lead to observable correlations. The persistent failure of this class 
of generators in the gap tests may lead to problems in some applications, e.g. in 
lattice simulations 



Finally, we would like to emphasize the importance of extensive testing such as pre- 
sented here before using any new pseudorandom number generator. Even a good 
algorithm can be corrupted by a poor implementation, as we have previously demon- 



strated |18| . Hence, a good amount of scepticism towards pseudorandom number 
generators without extensive test results seems prudent. It should also be noted that 
even when no statistical or bit level correlations are found, direct physical tests of 
random number generators should be used to reveal possible "hidden" correlations 

SB- 
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Table captions 



Table 1. List of the statistical tests. Numbers refer to the choice of parameters in 
Ref. 



Table 2. Results of the statistical tests. II, 12 and 13 refer to RCARRY, its 
corrected version, and ADCARRY, respectively. Results for RCARRY are from 
Ref. |TB). Depicted numbers are the observed descriptive levels 5 + and 5~ of the 



test variables K + and K~, respectively. A generator was considered to fail the test if 
the descriptive level was less than 0.05 or more than 0.95. Single, double and triple 
consequtive failures are indicated by single, double, and bold lines, respectively. The 
numbers shown are from the first run only. 

Table 3. Results of the bit level (i-tuple test. The bits marked failed have failed 
the test twice. See text for details. 



Figure captions 

Figure 1. 24 bit binary representation of random numbers produced by implemen- 
tation RCARRY (II) on a 120 x 120 matrix. 

Figure 2. 24 bit binary representation of random numbers produced by implemen- 
tation 12 of RCARRY on a 120 x 120 matrix. 

Figure 3. 24 bit binary representation of random numbers produced by ADCARRY 
on a 120 x 120 matrix. 
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Table I 



lest 


Test method 


number 




1 


Equidistribution test (1) 


2 


Equidistribution test (2) 


3 


Serial test in 2 dimensions 


4 


Serial test in 3 dimensions 


5 


Serial test in 4 dimensions 


6 


Gap test (1) 


7 


Gap test (2) 


8 


Gap test (3) 


9 


Maximum of t test, t — 5 


10 


Maximum of t test, t = 3 


11 


Collision test (1) 


12 


Collision test (2) 


13 


Collision test (3) 


14 


Runs-up test 
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Table III 



Implementation 


Failing 


Number of 




bits 


"good" bits 


11 


25-31 


24 


12 


25-31 


24 


13 


23-31 


22 
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